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In two-dimensional quantum site-percolation square lattice models, the von Neumann entropy is 
extensively studied numerically. At a certain eigenenergy, the localization-delocalization transition 
is reflected by the derivative of von Neumann entropy which is maximal at the quantum percolation 
threshold p q . The phase diagram of localization-delocalization transitions is deduced in the extrap- 
olation to infinite system sizes. The non-monotonic eigenenergies dependence of p q and the lowest 
value p q ~ 0.665 are found. At localized-delocalized transition points, the finite scaling analysis for 
the von Neumann entropy is performed and it is found the critical exponents v not to be universal. 
These studies provide a new evidence that the existence of a quantum percolation threshold p q < 1 
in the two-dimensional quantum percolation problem. 

PACS numbers: 71.30.+h, 03.67.-a, 72.15.Rn 
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I. INTRODUCTION 

The Anderson models and the quantum percolation 
(QP) model 3 -! 4 - are two important theoretical models that 
are used to study electron localization properties in dis- 
ordered systems. In the two models, the on-site energy 
randomness and the geometric randomness are consid- 
ered, respectively. For the Anderson model, there are ex- 
tensive studies and definitive results 2 -, while for the QP 
model, there are many open issues even today. 

The main concern in QP problems is to locate the QP 
threshold p q (accessible site concentrations by quantum 
particles) below which the electron is localized with prob- 
ability one. For the Anderson model and the QP model, 
there is a consensus on the existence of localization- 
delocalization transitions (LDTs) in three dimensions 
(3D)£i&££. For Anderson models, according to the one- 
parameter scaling theory 9 , LDTs do not occur at and be- 
low two dimensions (2D). However, whether the scaling 
theory is suitable for the QP model is an open question 1 ^ 
and even whether LDTs exist in 2D QP models at p q < 1 
is less clear—. Studies such as the scaling work based on 
numerical calculations of the conductance^ and transfer 
matrix methods with finite-size scaling analysis^ 3 - showed 
no evidence for LDTs. At the same time, there are 
many studies that claim LDTs exis t 14 ' 15 ' 16 ' 17 ' 18 . How- 
ever, the values of QP threshold p q obtained by different 
methods are not consistent. For example, for 2D quan- 
tum site-percolation square lattice models, Odagaki, et 
ai, obtained p q ~ 0.59 by a Green's function method^ 4 -. 
Koslowski and von Niessen gave p q ~ 0.70 with the Thou- 
less number based on the Thouless-Edwards-Licciardcllo 
methodic. Srivastava and Chaturvedi showed p q ~ 0.73 
with the participation used by the method with equations 
of motion^. Odagaki and Chang found p q ~ 0.87 using 
a real space renormalization-group method 1 ^. Ragha- 
van obtained p q ~ 0.95 by mapping a 2D system into 
a one-dimensional system 1 ^. Very recently, Islam and 



Nakanishi suggested that p q depending on particle en- 
ergies by calculating the transmission coefficient for 2D 
bond-percolation models 1 -^. 

In the mean time, metal states or metal-insulator tran- 
sitions are observed experimently in a variety of dilute 
two-dimensional electron and hole systems 1 ^. At the 
same time, the unusual transport properties of novel 
materials 1 ^, such as metal-insulator transitions happen 
in perovskite manganite films^S, and in granular metals 21 -, 
and minimal conductivity in undoped graphenes^ 2 -, may 
be explained by 2D QP models. Therefore these give 
additional motivations to study LDTs in 2D QP models. 

On the other hand, quantum entanglement, which 
attracting much attention in quantum information 2 ^, 
has been extensively applied in condensed matter 
physics^^^i^li^i^i^ ^^^^. For example, quantum 
entanglement measured by the von Neumann entropy 
was studied in the extended Hubbard mode l 26 ' 27 , in 
quantum small-world networks 2 ^, in two interacting par- 
ticle systems 3 ^, in the extended Harper mode l 27 ' 31 , in 
three dimensional Anderson models 3 - 2 -, in the integer 
quantum Hall system 3 - 2 -, and in spin model s 33 ' 34 . It is 
found that the von Neumann entropy shows singular be- 
haviors at quantum critical points(QCPs). The deriva- 
tive of von Neumann entropy has very good finite size 
scaling behaviours close to QCPs even for quite small sys- 
tem size o 32 ' 33 ' 34 . Therefore it becomes a powerful method 
to quantify QCPs in various systems. 

In this paper, with the help of the von Neumann en- 
tropy we present a detailed numerical study of LDTs in 
the 2D quantum site-percolation model . Our studies 
show that a quantum site-percolation threshold p q < 1 
exists in the 2D QP problem. In the next section the 
QP model and the definition of von Neumann entropy 
are introduced. In Sec. IHII the numerical results are pre- 
sented. And we present our conclusions and discussions 
in Section IIVI 
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II. QUANTUM SITE-PERCOLATION MODEL 
AND VON NEUMANN ENTROPY 

A. Quantum site-percolation model 

Let us consider one-electron tight-binding Hamiltonian 
with diagonal disorder defined on square lattices of sites^ 
The on-site potential £ can be drawn from the bimodal 
distribution 

P&) = pSfo ~ U) + (1 - P)5& - (1) 

In the limit £b ~ £a — ► °Oj the electron moves only on 
a random assembly of A-lattice points. Without loss of 
generality we choose £4 = 0. In the situation the A- 
site occupation probability is p and the corresponding 
quantum site-percolation Hamiltonian reads 

Haa = -t (4* + H - c -)> ( 2 ) 

(ij)eA 

where the summation extends over nearest-neighbor A- 
sites only. 

B. von Neumann entropy 

The general definition of entanglement is based on the 
von Neumann entropy^!. The generic eigenstate \a) for 
Hamiltonian with eigenenergy e a is the superposition 

|a)=E^H) = E^ c llO), (3) 

ieA i£A 

where |0) is the vacuum and ipf is the amplitude of wave 
function at ith site. For an electron in the system, there 
are two local states at each site, i.e., and |0}^, cor- 
responding to the state with (without) an electron at 

the zth site, respectively. The local density matrix pi is 
defined^^^^^ by 

ft = «4|l>«<l| + (l-«i)|0> ii <0| ) (4) 

where Zj = (a|cjcj|a) = \ipf 2 is the local occupation 
number at ith site. Consequently, the corresponding von 
Neumann entropy related to the ith site is 

Ki = ~ Z i l0 §2 z i ~ (! - z lo g2(! - z i)- ( 5 ) 

For nonuniform systems, the value of E^ depends on 
the site position i. At an eigenstate a, we define a site- 
averaged von Neumann entropy 

1 N 

E *=nY. E ^ ( 6 ) 

i=i 

were N is the number of A-lattice sites. The definition 
© shows that for an extended state that ib? — -4= for all 

i, E« = log 2 jf - (1 - log 2 (l - w i log 2 N at 



N — ► oo, and for a localized state that ipf = Sa ( i$ is a 
given site ) , = 0. In the present paper all the values 
of E" is scaled by jj log 2 N. From the two examples, 
we know the scaled E™ is near 1 when eigenstates are 
extended, and near zero when eigenstates are localized. 
Henceforth, we omit "scaled" for simplicity. 

For a random system the site-averaged von Neumann 
entropy E" should be further averaged over different re- 
alizations of disorder. The resulting quantity, the disor- 
der averaged von Neumann entropy denoted by (E") for 
eigenenergy e a , which is defined as 



N N 

i=l i=l 

where X is denoted as random average, K is the number 
of disorder realizations. In practice, (E") is the average 
values over a small window A around an energy value 
s, i.e., e a e [s - A/2,e + A/2]. We ensure that A is 
sufficiently small, and at the same time there are enough 
states in the interval A. Here A = 0.04 is chosen and 
other A give similar results. 

Another quantitative measure that is widely used to 
characterize localization is the participation ratio(PR)^, 

N 

defined by £ Q = (N l^f | 4 ) -1 ; which gives the ratio of 

lattice sites occupied by particles to all lattice sites at an 
eigenstate a. For the above extended and fully localized 
states, £ a = 1 and 1/N, respectively. Generally speak- 
ing, the larger £ a is, the more delocalized the eigenstate 
is. Similarly, the disorder averaged PR (£ Q ) is also av- 
eraged over the same small energy window. Henceforth, 
the disorder averaged von Neumann entropy (E£) and 
the disorder averaged PR (£ a ) are simplified to the von 
Neumann entropy (E£) and the PR (£ Q ), respectively. 



III. NUMERICAL RESULTS 

In numerical calculations, we directly diagonalize the 
eigenvalue Eq.([2]) with the periodic boundary condition 
and obtain all eigenvalues e a and the corresponding 
eigenstates \a). Without loss of generality, the hopping 
integral t is taken as units of energy. From formulas ([5][7]), 
we can obtain the von Neumann entropy (E"). We con- 
sider systems of linear size L = 20, 30,. ..,60 (measured 
in units of the lattice constant) and having N ~ L 2 A- 
lattice sites in total. The corresponding number of dis- 
order realizations K = 2500, 2000, 500, respectively. 
More realizations give similar results. 

Before discussing possible localized-delocalized transi- 
tions let us investigate the behavior of the von Neumann 
entropy (E%) at different site occupation probability p. 
The {E") as functions of eigenenergies e a are plotted in 
FigsQTa)-(c) a t L = 20,40 and 60 for p = 0.2,0.7 and 
0.9, respectively. It shows that at the same p, the values 
of {E") depend on eigenenergies e a and system sizes L. 
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FIG. 1: The von Neumann entropy (E") varying with 
eigenenergies e a at L — 20, 40 and 60 for site occupation 
probability (a)p = 0.2,(b)p = 0.7 and (c)p = 0.9, respectively. 
(d)The corresponding PR (£ a ) at p = 0.2,0.7, and 0.9 for 
L = 60. 
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FIG. 2: The relation between the {E") and the corresponding 
PR {^a) at L = 60 for states shown in Figfj] 



For p = 0.7 and 0.9, the (E") are relatively small for 
eigenstates near the band edges, while relatively large 
for eigenstates near the band center except eigenstates 
very near the band edges. Comparing the values of (E") 
for the three p at the same L, on the whole, all (£"") 
are relatively small for p = 0.2 and relatively large for 
p = 0.9. Compared to our results, the PR are also plot- 
ted in FigQJd) for L = 60 at the three p. The variations 
of the PR (£ a ) with respect to e a are similar to that of 
(E") in FigsQJa)-(c). At the same time, the {E") versus 
the corresponding £ Q are plotted in Figf2l It shows that 
the (E") increases monotonously with the so the 

von Neumann entropy can reflect the localization prop- 
erties of eigenstates in the QP model. 

To exact locate localized-delocalized transition points 
and the corresponding critical exponents with a finite- 
size scaling analysis, in the follows, we study the von 
Neumann entropy (E") changes with the site occupa- 
tion probability p for different eigenenergies . We found 
that the general trend for all von Neumann entropy (E£) 
curves at different eigenenergies are similar except the 
one at the band center , so we will discuss only two ener- 
gies as examples, one that is away from the band center 
and one very close to the band center. 



A. Eigenenergies away from the band center 

We first study the von Neumann entropy (E") and re- 
lated quantities for eigenstates with eigenenergies e a £ 
[-0.84 - 0.02, -0.84 + 0.02]. The literatures^iii^^iii 
agree on the QP threshold p q > p c — 0.593 if p q exits, 
where p c is the classical percolation threshold^. There- 
fore, we study the QP model at the site occupation prob- 
ability p beginning from 0.4, which is far smaller than the 
lower bound of p q . 

In FigPJJa), we show the dependence of (E") on 
the site occupation probability p at system sizes L = 
20,30, ...,60, respectively. It shows that (E") monotoni- 
cally increases as p becomes larger. For a certain system 
size, when p is small, e.g., p = 0.4, the eigenstates are 
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localized and (E") is small. When p = 1.0, the model 
shown in Eq.([2]) is a two-dimensional periodic potential 
system. Due to the Bloch theorem the eigenstate of 
a tight-binding electron on a local regular lattice is al- 
ways in the extended state. At the situation, (E™) is 
largest. All these reflect the trivial derealization effect 
of p, which is similar as that studied in quantum small- 
world network models 2 ^.. All data shown in Fig[3]Ja) are 
well fitted with nonlinear Boltzmann functions for vari- 
ous system sizes. According to the fitted lines, we plot 
the derivative d{E")/dp varying with p in Fig[3Jb). ft 
shows there is a peak in the derivative at a certain p, 
which is denoted by p max - The maximal derivative and 
Pmax increase with the system sizes L, respectively. It 
is believed that the von Neumann entropy may be non- 
analytic at a quantum phase transitions and can reflect 
various quantum critical point o 26 ! 32 . Therefore, that the 
derivative is maximal at a certain position p max , can be 
as a signature of LDTs of electron state a 27 i 29 . 

To study the LDTs at the QP threshold p q , one needs 
to investigate the behavior of systems in the thermody- 
namic limit. However, in most cases this is not possible 
in numerical methods 3 - 7 -, and therefore, similarly as in 
Refl27l an extrapolation method is chosen. Fig(3]Jc) shows 
the scaling behavior of the p max - The p q in the thermo- 
dynamic limit can be obtained by 1/L — > and Figj3jc) 
shows p q ~ 0.676 at the situation. We denote the deriva- 
tive of von Neumann entropy at Vmax as d(E") / dp\ max 
. Following the Refs. HH and [3J, the finite size scal- 
ing is performed for the function 1 — exp(d(E") /dp — 
d(E%) I dp\ max ) with respect to L 1 /" {p-p max ). The result 
is presented in Fig(3Jd). It shows numerical results ob- 
tained from various system sizes approximately collapse 
on a single curve with the critical exponent v ~ 2.52. 



B. Eigenenergies near the band center 

In the following, we discuss the von Neumann entropy 
(E") for eigenstates with eigenenergies e a E [0 — 0.02, + 
0.02]. The von Neumann entropy {E") and the corre- 
sponding derivative d(E")/dp varying with the site oc- 
cupation probability p are shown in Figs [Ha) and (b), re- 
spectively. Figdja) shows that (E") first increases with 
p until to a plateau at p 6 [0.6,0.85], then continues 
to increase, which is quitely different from that shown 
in Fig[3Ja) for eigenenergies away from the band cen- 
ter. FigfJJb) shows that the derivative d(E")/dp dras- 
tically decreases near p ~ 0.6 and drastically increases 
near p ~ 0.85. 

With a real-space renormalization-group method, Oda- 
gaki and Chang observed three regimes of the electronic 
properties in QP models, which divided by the classical 
percolation threshold p c and the quantum percolation 
threshold p q {p q > p c )r^- When p < p c , electrons can- 
not tunnel between different isolated clusters and elec- 
trons are considered to be localized in the classical sense 
even in the quantum case. When p c < p < p q , due 
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FIG. 3: Some quantities for eigenstates with eigenenergies 
e a e [-0.84 - 0.02, -0.84 + 0.02]. (a) The von Neumann en- 
tropy {E") varying with the site occupation probability p and 
the lines are Boltzmann fitting, (b) The derivative d{E")/dp 
varying with p. (c) The line corresponds to the expected be- 
havior of p m ax for 1/L — > according to a second-order 
polynomial fitting, (d) The finite-size scaling analysis. The 
system sizes L = 20, 30, 60 and the arrow direction in (a) 
and (b) denotes the increasing of L. 
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(b) P 

FIG. 4: Some quantities for eigenstates with eigenenergies 
e a 6 [0 - 0.02, + 0.02]. (a) The von Neumann entropy (JS«) 
and (b) the d{E")/dp varying with the site occupation prob- 
ability p. The system sizes L — 20,30, ...,60 and the arrow 
direction in (a) and (b) denotes the increasing of L. The inset 
in (b) is for L = 60. 



to quantum interference effects, electrons cannot spread 
infinitely even there is an infinitely extended channel. 
The regime is called quantum localization regimes. When 
p > p q , electrons can spread infinitely and electron states 
are extended. They found p c = 0.618 and p q = 0.867. It 
is interesting that the von Neumann entropy (E£) dras- 
tically changes near the two p values. 



C. Phase diagram 

We have studied and extensively analyzed the von Neu- 
mann entropy for all the other eigenstates. As the system 
in Eq.([2]) is bipartite, e a and — e a are both eigenvalues 
of Haa—- Therefore we will restrict our investigation 
to one (left) half side of the band, i.e.,— 4 < e a < 0. 
FigEIa) presents the phase diagram of LDTs in the e a —p 
plane. The values of p at the LDTs are the QP thresh- 
old p q , which is obtained by the extrapolation method 
similarly as that shown in Fig[3Jc). The trend for p q 
varying with e a is similar as that for 2D quantum bond- 
percolation models on square lattices^ and 3D quantum 
site-percolation models on simple cubic lattice^. In de- 
tail, we observe the nonmonotonic dependence of the val- 
ues of p q on eigenenergies e a . Notice that in the region 



— 2t < e a < — 0.5i the QP threshold p q are nearly con- 
stant and a weak maximum at e a ~ —t, which are very 
similar as that found in Ref.[f|. That a weak maximum 
for p q at e a ~ —t has also been found and discussed 
in Ref.[8[, which may be due to the existence of von 
Hove singularity at the energy^. We find all the QP 
threshold p q are greater than the classical percolation 
threshold p c ~ 0.593 and the lowest value p q ~ 0.665. 
As shown in Figj5ja), the phase diagram is very consis- 
tent with the mobility edge trajectory shown in Ref.[l5T| 
for the same model obtained by the Thouless-Edwards- 
Licciardello method. 

FigOJb) presents all the critical exponents v versus 
eigenenergies e a according to the finite scaling analy- 
sis similarly as that shown in Fig(3Jd). We find the 
values of v depend on eigenenergies e a . In detail, in 
the region — 2t < e a < — 0.5i where p q are almost 
constant (see Figj5ja)), most of the values of v are dis- 
tributed in a relatively narrow interval [2.6,3.2]. From 
the band edge to e a ~ — 2t, v increase with e a , while 
from e a ~ — 0.5t to the band center, v decrease with 
e a . All v arc larger than 2/D(D = 2), which satisfies 
the assumption that v must satisfy the bound v > 2/D 
for random systems^. Though v has been extensively 
studied in 3D QP models&iii, to our best knowledge, 
there are few works to study v in 2D QP models except 
in Ref. [l7l where the critical exponent for correlation 
lengths v ~ 3.35. The varying of the critical exponent v 
with the QP threshold p q is plotted in FigOJc) . It shows 
near the band edge and near the band center the relation 
between v and p q is linear, respectively, but the linear 
relations are different. This may be caused by the dif- 
ferent symmetry and/or degeneration at the two energy 
regions^. 



IV. CONCLUSIONS AND DISCUSSIONS 

In this paper, we have detailed studied the von Neu- 
mann entropy (E") varying with eigenenergies e a and 
accessible site concentrations p by quantum particles 
in two-dimensional quantum site-percolation models on 
square lattices. 

For Eigenenergies away from the band center, we de- 
termine the QP threshold p q by the derivative of von 
Neumann entropy is maximal at the point. Based on 
this, we give the phase diagram of LDTs in the e a — p 
plane, which is very consistent with the mobility edge 
trajectory shown in Ref. [l5j . From the phase diagram, 
we observe the non-monotonic eigenenergies dependence 
of p q and the lowest value p q ~ 0.665. At the same time, 
the finite-size scaling analysis is performed at all LDTs 
points. To the best of our knowledge, it is the first time 
to obtain all the critical exponents v at the whole energy 
space for 2D QP models. We find the critical exponents 
v depend on eigenenergies e a . 

For eigenenergies near the band center, the variations 
of the von Neumann entropy (E£) with respect to p 
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All our numerical results show that there is p q < 1 in 
2D QP models and p q depends on eigenenergies. The de- 
bates on the values of p q may be partially due to different 
energies treated in literatures. 
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